Thermal conductivity in a mixed state of a superconductor at low magnetic fields 
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We evaluate accurate low-field/low-temperature asymptotics of the thermal conductivity perpen- 
dicular to magnetic field for one-band and two-band s-wave superconductors using Keldysh-Usadel 
formalism. We show that heat transport in this regime is limited by tunneling of quasiparticles 
between adjacent vortices across a number of local points and therefore widely-used approximation 
of averaging over the circular unit cell is not valid. In the single-band case, we obtain parameter-free 
analytical solution which provides theoretical lower limit for heat transport in the mixed state. In 
the two-band case, we show that the heat transport is controlled by the ratio of gaps and diffu- 
sion constants in different bands. Presence of a weaker second band strongly enhances the thermal 
conductivity at low fields. 



I. INTRODUCTION 

The thermal conductivity of a metal in the super- 
conducting state is markedly different from its value 
in the normal state. The physical reason is that the 
Cooper pairs give no contribution to the transport of 
heat. Therefore heat transport occurs solely due to the 
quasiparticle excitations controlled by the energy gap. 
Due to this reason, measurements of the thermal con- 
ductivity were extensively used during last decades as a 
tool to probe the energy gap symmetry in various new 
superconducting materials. 

Application of a magnetic field H provides a way to 
generate quasiparticles in a type-II superconductor. In a 
s-wave superconductor at T <C T c , the quasiparticles are 
localized near the vortex cores and thermal conduction 
perpendicular to the magnetic field is due to tunneling 
between adjacent vortices. Thermal conductivity in the 
mixed state of superconductors was extensively studied 
in the past experimentally^— and theoretically 4 ^. Early 
theoretical work 4,5 addressed superconductors in diffu- 
sive regime (dirty limit) in the mixed state with a tem- 
perature gradient applied transverse to the magnetic flux. 
The electronic thermal conductivity was calculated using 
the quasiclassical time-dependent superconductivity the- 
ory, assuming homogeneous temperature gradient, using 
the circular-cell approximation, and averaging over the 
unit cell. This approach was later extended to study 
heat transport in various unconventional superconduc- 
tors. Theory for d-wave superconductors was developed 
in Ref. for zero magnetic field and in Refs. and 
8 for the mixed state. More recently, with the dis- 
covery of magnesium diboride (MgB2), theory was ex- 
tended to the case of multiband superconductivity^. Dis- 
covery of pnictides, with possibly unconventional pair- 
ing mechanism due to spin fluctuations leading to the 
so-called s± pairing state, motivated extension of heat 
transport theory^ taking into account multiband su- 
perconductivity and resonant interband impurity scat- 
tering. Many recent experimental studies reconsider old 
superconductors, such as NbSe? 1 ^, and addressed novel 



superconductors, such as borocarbides (LuNi 2 B 2 C) 11 , 
Sr 2 RuO^, CeYbH, heavy- fermion compounds Celrln 5 , 
CeCoIn 5 and UPt 3 13 , MgB 2 15 , pnictides 1 ^ and iron- 
silicides 18 . Being very sensitive to the gap structure, 
thermal transport at low magnetic fields varies consid- 
erably for different compounds. 



In this paper we reconsider more accurately the prob- 
lem of the thermal conductivity across the magnetic fields 
for s-wave superconductors. The case of s-wave supercon- 
ductor is a standard reference, all other situations are 
compared with this case. Surprisingly, an accurate result 
for the low-field asymptotics of the thermal conductiv- 
ity was never derived. The widely-used circular unit cell 
approximation does not give correct result in low mag- 
netic fields because it misses essential physics. Namely, 
as the local thermal conductivity is strongly inhomoge- 
neous, the heat transport is limited by tunneling between 
adjacent vortices across certain local points in the vortex 
lattice unit cell (bottlenecks). This leads to general low- 
field asymptotics of the electronic thermal conductivity, 
k oc exp(— P^f B C 2/ B), where -B c2 is the upper critical 
fields Surprisingly, the theoretical value of the numeri- 
cal constant j3 is not available neither for clean nor for 
dirty s-wave superconductors. For clean case, we provide 
estimate for this numerical constant using asymptotics of 
the Bogolyubov wave functions of the localized states at 
zero energy and microscopic value of the upper critical 
field. In the dirty case we were able to perform more 
quantitative analysis using the Keldysh-Usadel formal- 
ism. We calculate the thermal conductance at low tem- 
perature and low magnetic field for single- and two-band 
superconductors in the dirty limit. In this regime we ob- 
tain parameter-free analytical solution in a single-band 
case which provides theoretical lower limit for heat trans- 
port in the mixed state. We find that in dirty case the 
low-field thermal conductivity is drastically suppressed 
in comparison with clean case. Further, we generalize 
the developed formalism to a two-band superconductor, 
taking MgB 2 as an example. 
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II. TUNNELING OF QUASIPARTICLES 
BETWEEN VORTEX CORES IN CLEAN 
ISOTROPIC SUPERCONDUCTOR 

The electronic transverse thermal conductivity in 
mixed state at low temperatures and fields is determined 
by the probability quasiparticle tunneling between the 
cores of neighboring vortices. In this section we eval- 
uate this quantity with exponential accuracy for clean 
isotropic superconductors. Even though this estimate is 
very straightforward and could be done long time ago, to 
our great surprise, we did not succeed to find it in the 
literature. 

The Bogolyubov wave function of the localized state 
in the vortex core at E — decays asi£ 



*(r) oc exp(-r/£ A ), £a = v F /A 



(1) 



where A is the superconducting gap and vf is the Fermi 
velocity. Therefore, the probability of tunneling between 



the vortex cores separated by distance a 
can be estimated as 



P oc | * (a) | 2 oc exp (- 2a/£ A ) = exp 




2$ /V3B 



(2) 



The upper critical field B C 2 for a clean isotropic super- 
conductor at T — is given by2£ 



B C 2 = 



e 2 $o 



VF 

2irT r 



(3) 



47 2tt$ 

with 7 = exp(0.5772) = 1.781 being the Euler constant. 
Using also the BCS relation ttT c — 7 A, we evaluate 



2$n 



' I6n 



1.486. 



This allows us to represent the tunneling probability (|5J) 



MS 



Pocexp 



16tt B 



c2 



B 



;exp -1.486 




(4) 



This result also determines the low-field asymptotics of 
the electronic thermal conductivity of clean isotropic su- 
perconductor. The constant in the exponent is obviously 
sensitive to anisotropy of Fermi surface. A more quanti- 
tative analysis which would include also evaluation of the 
preexponential factor requires much more complicated 
microscopic kinetic theory for clean limit. In the follow- 
ing sections we make quantitative calculations of thermal 
transport at low fields for dirty superconductor. 



III. THE FORMALISM: THERMAL 
TRANSPORT IN DIRTY SUPERCONDUCTORS 
WITHIN QUASICLASSICAL KELDYSH-USADEL 
MODEL 

Our study is based on the quasiclassical Keldysh- 
Usadel formalism 2 ^— which was developed to describe 



nonequilibrium properties of dirty superconductors. Be- 
low we will reproduce the main relations of this formal- 
ism needed for our derivations. Within this formalism a 
superconductor is described by the Green's function 



G = 



G R G K 
G A 



(5) 



where in dirty limit the retarded (advanced) Green's 
functions G R ^ satisfy the Usadel equation 24 



hDV(G R WVG R W) = \iEr 3 + A,G fl ( A ) 



(6) 



Here E is quasiparticle energy, D is the electronic diffu- 
sion coefficient, A is the pair potential 



A = 



A 

A* 



(7) 



qK = qRj _ jqA^ j ig the distribution function, 
G R = t 3 G + nF, G A = -f 3 G fl tf 3 , where G and F are 
the normal and anomalous Green's functions and fi are 
Pauli matrices. 

Thermal current is given by 



NnD 



>th 



ETt 



t 3 (g r VG k + G k \7G a ) dE, (8) 



where No is the normal density of states the Fermi level. 
Writing f = f L i + f T T 3 , where fa and fr are odd and 
even in energy components of the distribution function, 
one can rewrite the thermal current in the form 2 ^ 



Jth = N, 



J E [D L (E) V/ l (E) + ImJ B / T ] dE, 



(9) 



where 



D T = D 
D L = D 



(ReG) 2 + (ReF) 2 
(ReG) 2 - (ImF) 2 



are the energy-dependent spectral diffusion coefficients 
and ImJ^; is the spectral supercurrent given by ImJ^ = 
iTr^G^VG^ - G A VG A )} = ImF R ReF R V X , where X 
is the superconducting phase. In the mixed state all 
quantities are spatially inhomogeneous (coordinate de- 
pendences are dropped for brevity). Functions /l and 
Jt satisfy the following kinetic equation o 22 i 23 

V(P T V/ T ) + ImJ B V/ L = 2i?/ T , (10) 
V (D L Vf L ) + ImJ B V/ T = 0. (11) 



where R 



In thermal equilibrium 



A(G R + G A ) 

Jl = tanh (i?/2fcsT) and fa = 0. The thermal conduc- 
tivity components are defined as 



K a = (Jt/i,a>/(V Q T), 



(12) 
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FIG. 1. Left: Graylevel map of the local thermal resistance 
oc 1/Dl(0, r) in the vortex lattice. Light regions correspond 
to low thermal resistance. Arrows illustrate heat flow in the 
bottleneck regions. Right: Three-dimensional plot of the local 
thermal resistance in the bottleneck region marked by rect- 
angle in the left picture. 



where (VT) is the average temperature gradient and a = 

(x,y,z). 

In the following, we shall use the standard 0- 
parametrization, G R (r) = T3cos[0(r)] + fi sin[0(r)] in 
which the Usadel equation has the following for m 26 ' 27 

L> ( V 2 - p 2 cos sin 0) + 2A cos + 2iE sin = (13) 

where p(r) = V<fi — (2-7r/<I>o)A is superconducting mo- 
mentum (within circular cell approximation p = l/r — 



r jr\ 
the form 



= y/Qo/irH). The selfconsistency equation has 



the boundary regions between the cells, see Fig. [TJ This 
means that thermal transport occurs via "bottlenecks", 
saddle points of Dl(E,y). Our purpose is to evaluate 
average thermal conductivity limited by these bottle- 
necks. Since in the vicinity of bottlenecks, the spectral 
supercurrent ImJs vanishes by symmetry, the expres- 
sion for local energy current simplifies and has the form: 
J th (r) ksNoJ dEED L \7f L . 

We start with evaluation of the diffusion constant 
Dl (E, r) which is determined by the real part of the 
Green's function 8(E,r). For an isolated vortex, at dis- 
tances r>( from its core we can present A and 9 as 
0(r) = O + 0(r), A(r) = A + A(r) where A and 
0o are the equilibrium values at zero magnetic field, 
tan 0o = iAo/E and 0(r) = 9 r (r) + %9i (r) is small cor- 
rection. For E <C A , 0o ~ 7r/2 + i_E/A , therefore the 
energy-diffusion constant in this region Dl(E, r) « D9 2 . 
The real and imaginary parts of obey the following 
equations 

D E l + A % 2 9 r - 2\lA 2 -E 2 9 r = 0, (15) 



D\I l 9 r 



DXJ 2 9 1 - D 



E* -A 2 



E 2 



E 2 



A 2 , 



E 2 9 t 



2EA 



E 2 



- D 



EAn 



A 2 



E 2 



P 



(16) 



where p = l/r is the gauge-invariant phase gradient. For 
an isolated vortex 9 r and 9i have qualitatively different 
behavior at large distances for E < Aq: 9% decays as 
p 2 cx l/r 2 and 9 r decays exponentially 



exp(— k^r) 



(17) 



j< j 1 / ^ 

A In — = — sin 



(14) 



The diffusion coefficients in the 0-parametrization 
are given by the expressions Dt{E,t) = 
D cosh 2 (Im[0(J5, r)]), D L (E, r) = D cos 2 (Re[9(E, r)]). 



IV. THERMAL TRANSPORT IN THE VORTEX 
STATE AT LOW FIELDS 

We consider a superconductor in low magnetic field, 
B -C B c2 . We assume an ideal triangular vortex lat- 
tice of vortex lines and study thermal transport in the 
direction perpendicular to the field. The local ther- 
mal conductivity is mostly determined by the diffusion 
constant Dl(E, r) at low energies. This quantity is 
very inhomogeneous in the vortex state. It has max- 
ima at the vortex cores and rapidly drops away from the 
cores reflecting localization of quasiparticles in the core 
regions. This means that the local thermal resistance 
cx 1/Dl(0,t) is maximal at the boundaries of the lattice 
unit cell. In such situation the temperature is mostly 
homogeneous within the unit cells and only changes in 



V2(A 2 Q ~E 2 ) 1/4 / 



D. 



Evaluation of numerical constant requires solution of full 
nonlinear problem. Numerical solution by the method 
described 27 provides C v w 4.2. In the vortex lattice away 
from the core regions 9 r can be represented as a sum of 
contributions from individual vortices 



R) 



=5> 



R 



R 



exp(— k(\r — R| 



Rl 



(18) 

where R are the vortex coordinates. In particular, for 
triangular lattice R = (a/2 + ma + na/2, n^/i/2) where 

a = \J 2$o/v / 3-B is the lattice constant and m and n are 
integers. 

At small field thermal transport is determined by en- 
ergy flow via bottlenecks, saddlepoints of Dl(E,t) at 
boundaries of the lattice unit cell. Near the bottleneck 
point, (a;, y) = (0, 0), we can keep only contribution from 
two neighboring vortices located at (±a/2, 0) which gives 

2C V 



'lat,- 



(r) 



Vha~j2 



exp (— k^a/2) cosh (k^x) exp(— k^y /a) 

(19) 
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and 



8C 



Dl(E, r) sa D- — - exp (—k^a) cosh (k^x) exp(—2k^y ft 



k^a 



(20) 



Using also quasiequilibrium approximation 
for the gradient of f L (E,r), V x f L {E,r) « 
-cosh~ 2 (E/2k B T)(E/2k B T 2 )V x T(r), we obtain the 
heat flow near the bottleneck Jth, x (r) ~ — K x (r)W x T(r), 
where the local thermal conductivity n x (r) is given by 



« x (r) « 6(%DNq 



dEE 



exp ( 



-fetal 



cosh 2 (k^x) exp 



2k^ 



E/2k B T 2 



cosri (E/2k B T) 



r 



At low temperatures T < Ao£/a the main contribution 
to the energy integral comes from the region E <T. This 
allows us to neglect the energy dependence of k% (E) and 
replace k^(E) -> fc^ = y / 2A Q /D = 1/£a- In this case, 
using x 2 cosh _2 x dx — 7r 2 /6, we obtain 



Kx(r) 



^Cgfc|^n^ eXP . ( ~ %0a) 



HO' 



cosh (fc^o^) exp I — 



2k (0 y< 



Near the bottleneck region the local thermal conductivity 
has a form n(x,y) = KoF x (x)F y (y), where the function 
F x {x) — cosh 2 (k^ox) has minimum at x — and F y (y) = 
exp (— 2k^ y 2 /a) has maximum at y = 0. The total flow 
per unit length along the field through the bottleneck is 
given by 



hh ~ koF x (x) 



F y {y)dy 



(22) 



vortex-lattice state 
k/T 



IOv^^C 2 2 



3-y/3 



10V27r 5 / 2 C 2 
3V^ 



exp (-fcg q) 
^/k^a 

exp (-J{8n/V3)B c2 /B 



klN D- 



((8ir/V3)B c2 /B 



1/4 



(26) 



(21) where we used relations B c2 = <f>o^|o/47r and fc{o« 



n/V3)B c2 /B. 

Introducing the thermal conductivity in the normal 
state K W = ^k%N DT, we rewrite the final result in 
the form 



K 
KN 



— w 10 



130 



exp— J(87r/V3)B c2 /B 



((8Tt/V3)B c2 /B 



1/4 



(27) 



Energy conservation requires that / t /j has to be x- 
independent. Therefore, the temperature drop AT across 
the bottleneck can be evaluated as 



AT: 



ST 00 F x \x)dx 

K oJZo F v(y) d y 



Hh 



(23) 



meaning that the total thermal conductance through this 
region K — I t h/AT can be evaluated as 



This parameter-free analytical result provides theoretical 
lower limit for the heat transport in the mixed state in an 
isotropic dirty s-wave superconductor at low field. The 
constant 3.81 in the exponent is significantly higher than 
the constant 1.486 which we evaluated for the clean case 
meaning that the scattering drastically suppresses the 
quasiparticle thermal conductivity at low fields. 



V. DISCUSSION OF EXPERIMENT 



K = k q 



IT 00 F y (y)dy 

JZ a F x 1 (x)dx 3V2 



n5/2C h%TN D eXp{ - k ^ a) 



Vafcfo 

(24) 

Evaluating the total average energy flow density 

J th = 7 " + W4 = *K**, (^) 
aV3/2 2^3 a ' K ' 

we obtain the final result for the low- field/low- 
temperature limit for the thermal conductivity in the 



Rather limited set of experimental data is available on 
electronic contribution to the thermal conductivity at low 
temperatures and low magnetic fields, since phonon con- 
tribution should be accurately subtracted in this regime. 
To our knowledge, the only experimental data on elec- 
tronic thermal conductivity which clearly demonstrate 
low-field/low-temperature behavior expected for an s- 
wave superconductor are rep orted for Nb in Refs. [l| and 

and for VsSi in Ref. |TJ In these experiments the 
samples were in the clean limit. In Ref. [H the elec- 
tronic thermal conductivity was fitted by the expression 



5 



K oc exp(— f3y/B C 2/B) for values of B up to about -B c2 /3, 
with (3 = 1.66. This number is slightly higher than our 
estimate of j3 = 1.486 in Eq. (|4|. This difference can be 
explained by the influence of impurity scattering. The 
shape of field dependence of thermal conductivity for 
VaSi reported in Ref. [l2|is in qualitative agreement with 
the predicted exponential dependence, however the quan- 
titative analysis was not made and the value of j3 was not 
explicitly extracted. In a field B = B C 2/20 the value of 
k « 2.5 • 10~ 3 kjv provided in Ref. [H exceeds our dirty- 
limit estimate, Eq. (f2"Tl) . by about two orders of mag- 
nitude. This discrepancy can be naturally attributed to 
the fact that measured VaSi samples were in the clean 
limit. The magnitude of thermal conductivity at not 
very small magnetic fields is in good qualitative agree- 
ment with calculations made for the clean limit using 
the Landau-level expansion and assuming homogeneous 
temperature gradient^. To our knowledge, there are no 
data available on electronic thermal conductivity of dirty 
s-wave superconductors in the low- field/low-temperature 
limit. In available measurements of alloys 1 and dirty Nb 
samples 3 - the thermal conductivity at low fields and tem- 
peratures is dominated by phonons and separating elec- 
tronic contribution is a challenging task. 



In the presence of two electronic bands with different 
energy gaps, low-temperature behavior of thermal trans- 
port is determined by the band with lower gap (7r-band) . 
Still, the generalization from single-band to two-band 
case involves not simply renormalization of the energy 
gap, but also correction to the asymptotic behavior of 
Green's function in the it -band and to the upper critical 
field. 

The Usadel equation for the 7r-band reads 

( V 2 ^ ~p 2 cos 6^ sin 9„) +2A W cos 6 n + 2iE sin 6 7T = 0, 

(28) 

where and A^r) are the diffusion constant and gap 
for the 7r-band. Similar to a single-band case, asymp- 
totics of the Green's function in the 7r-band at large dis- 
tance from the vortex core is given by 



^7T,r(^*) C*7r 



exp(— k^r) 



(29) 



with 



2\l/4 



The main exponential dependence of the zero-energy 
Green's function at the bottleneck point is oc exp(— k^a), 
where /c^q — &tt(0) = y / 2A 7I o/ D^, which gives 



VI. EXTENSION TO A TWO-BAND 
SUPERCONDUCTOR 



Here we extend the above formalism to a two-band 
superconductor. At present, the most established exam- 
ple of such system is MgB2, which is characterized by 
two electronic bands: 7r-band and cr-band, see, e. g. re- 
cent review Ref. [29|. The quasi-two-dimensional cr-band 
is characterized by stronger superconductivity than the 
three-dimensional tt -band. Heat transport in a two-band 
superconductor was studied theoretically in Ref. [§| assum- 
ing clean limit conditions and using the averaging over 
unit cell method^. In the dirty limit, theory of density 
of states in the mixed state and the upper critical field 
for MgB2 was developed in Refs. and HH. Below we 
extend the calculations of the heat transport presented 
above, to the case of a diffusive two-band superconduc- 
tor, taking MgB2 as an example. 



k^a = \J Bk/B, 

2A.rO 2$p 
At V3' 



(30) 
(31) 



Therefore, the magnetic field dependence of thermal con- 
ductivity is determined by the field scale B^ and can be 
presented in the form similar to Eq. (|27l) . 



KttN 



10 



2- a2 exp 



(-^BJB) 



{Bk/B 



1/4 



(32) 



where k w n = ^-k^N^D^T is the partial 7r-band contri- 
bution to the normal-state thermal conductivity. 

To proceed further, we have to find relation between 
the 7r-band field scale B n and the upper critical field B c2 
for a two-band superconductor. The upper critical field 
at low temperatures, T <C T c , is given by 31 



B c2 (0) = a c2 B^. 2 (0), a c2 = exp 



W 1 +W 2 -lnr :i 



'(Wi+Wa-lnr a 



Wi In r x 



(33) 



where indices 1 and 2 correspond to the a and 7r bands, A a p is the coupling-constant matrix, r x — T> a /T>^ 1 T> a 

are the diffusion constants in a and 7r bands, 



T(An - A 22 )/2 + V(An - A 22 ) 2 /4 + Ai 2 A 2 i s , . _ A B cs$o 

AiiA 22 - A 12 A 2 i ' c2U ~ 2irV a 



6 



is the single-band upper critical field for the c-band, and 
A BCS = ne'^ E T c w 1.764T C (j E w 0.5772 is the Euler 
constant). 

To make estimates for MgB2, we use the follow- 
ing coupling matrix element a 30 ' 31 : An s» 0.81, A22 ~ 
0.278, Ai2« 0.115, A 2i « 0.091, which gives Wi 0.088 
and W2 ~ 2.56. With such coupling matrix the two- 
band BCS model gives 3 ^ A^o ~ 0.3A ct0 ^ 0.177ttT c . 
Since the parameter Wi is small, typically the inequal- 
ity Wi I Lara I "C {W 2 — lnr x ) /4 is valid. In this case 
one can expand Eq. (|3"3"|) with respect to Wi and obtain 
simple result 



a c2 ~ 1 



Wi In r a 
W 2 - In r. 



(35) 



meaning that the upper critical field is close to B s c2 (Q) 
and is mostly determined by the coherence length of the 
cr-band. Using Eqs. (|3"Tj) , (j3"3")l , and (fM]) . we obtain the 
relation between and £? C 2(0) 



Bit 



8tt A x0 £> ff 
-\/3a C 2 Abcs T^-k 



Sca(0), 



(36) 



which presents the main result of this section. This scale 
has to be compared with the scale (8Tr/y/3)B c2 (0) for the 
single-band case. 

In order to determine the pre-exponential factor CV 
in Eqs. (|29|) and (|32|). we have to calculate the Green's 
function niT {r), which requires solution of the full two- 
band Usadel problem, as described in Ref. [3(1 An im- 
portant parameter is the ratio of diffusion coefficients 
in two bands r x . For illustration, we consider two 
cases here: r x = 1 and 0.2, for which the ratios of 
the coherence lengths in the two bands are ^/Co- = 
v /(D ;r /D CT )(A c7 o/A 7r o) = 1-83 and 4.1. For these two 
cases we compute B T sa 0.32(87r/v3)-B C 2(0), C v ~ 3.6 
for r x = 1 and B„ w 0.065(87r/VS)B c2 (0), CV « 2.9 for 
r x = 0.2. This gives 



and important fingerprint of the two-band superconduc- 
tivity in MgB 2 . This small scale was experimentally ob- 
served not only in the thermal conductivity^, but also 
in the specific heat 3 ^ and flux-flow resistivity^ 3 -. For the 
magnetic field applied along c-axis it is 3-5 times smaller 
than B c2 . In this case, the low-field regime described by 
Eq. d32) is expected at fields < B c2 /30 w 100G. Unfortu- 
nately, most experimental data of Ref. [ID are presented 
for higher magnetic fields. 

We shall note that in available MgB 2 single crystals 
estimates suggest that the a band is in the clean limit. 
However, our results should be qualitatively applicable 
even in this case. The reason is that dominant contribu- 
tion comes from the 7r-band which is in the dirty limit, 
as argued in Ref. [3^ where the low-energy DoS in the 
vortex state of MgB2 was calculated. 

In summary, we have readdressed the problem of the 
heat transport of a superconductor in the mixed state 
at low temperatures and low magnetic fields, going be- 
yond the circular unit cell approximation. In the clean 
limit we estimated the numerical constant /3 in the low- 
field asymptotics of the electronic thermal conductiv- 
ity, k oc exp(— Py/ B c2 /B), using the Bogolyubov wave 
functions of the localized states at zero energy. In the 
dirty limit we have performed quantitative analysis of 
heat transport using Keldysh-Usadel formalism and have 
shown that heat transport is limited by tunneling be- 
tween adjacent vortices across certain local points (bot- 
tlenecks) . In the isotropic s-wave superconductor we have 
obtained parameter-free analytical solution which pro- 
vides theoretical lower limit for heat transport in the 
mixed state. Based on this solution, one can conclude 
that low- field/low-temperature thermal conductivity in 
the mixed state is drastically suppressed by impurity 
scattering. We have extended our results to the case of a 
two-band superconductor, taking MgB 2 as an example. 
In this case, we predict an enhancement of heat transport 
with strong dependence on the ratio of gaps and diffusion 
constants in different bands. 



K I B 

1 130 



1/4 



:130 



B 



B c2 (0) 



1/4 



exp -0.98 



B c2 {0) 



B 



, r x = 0.2 



We can see that the field scale in the two-band case is 
strongly reduced in comparison with the single-band case 
leading to large enhancement of the thermal conductiv- 
ity at low fields. This reduction is mostly caused by the 
smaller energy gap in the 7r band. Another factor which 
may contribute is possible large value of the diffusion 
constant . The smaller field scale caused by the larger 
coherence length in the 7r-band is an established feature 



ACKNOWLEDGMENTS 

The authors would like to acknowledge useful dis- 
cussions with N.B.Kopnin, V.M.Vinokur, and I.I.Mazin. 
A.E.K. is supported by UChicago Argonne, LLC, opera- 
tor of Argonne National Laboratory, a U.S. Department 
of Energy Office of Science laboratory, operated under 
contract No. DE-AC02-06CH11357. This work also was 
supported by the "Center for Emergent Superconduc- 
tivity" , an Energy Frontier Research Center funded by 
the U.S. Department of Energy, Office of Science, Of- 
fice of Basic Energy Sciences under Award Number DE- 
AC0298CH1088. 



1 J. Lowell, J. B. Sousa, J. Low Temp. Phys., 3, 65 (1970). 2 W. F. Vinen, E. M. Forgan, C. E. Cough, and M. J. Hood, 



7 



Physica 55, 94 (1971). 

3 P. H. Kes, J. P. M. van der Veeken, and D. de Kierk, Journ. 
of Low Temp. Phys., 18, (1975). 

4 W. Pesch, R. Watts-Tobin, and L. Kramer Z. Physik 269, 
253 (1974). 

5 S. Imai and R.J. Watts-Tobin, Journ. of Low Temp. Phys., 
26 967 (1977). 

6 M. J. Graf, S-K. Yip, J. A. Sauls, and D. Rainer, Phys. 
Rev. B 53, 15147 (1996). 

7 A. B. Vorontsov and I. Vekhter, Phys. Rev. B 75, 224502 
(2007). 

8 C. Kubert and P. J. Hirschfeld, Phys. Rev. Lett. 80, 
4963(1998). 

9 H. Kusunose, T. M. Rice, and M. Sigrist, Phys. Rev. B, 
66, 214503 (2002). 

10 V. Mishra, A. Vorontsov, P. J. Hirschfeld, and I. Vekhter, 
Phys. Rev. B 80, 224525 (2009) 

11 E. Boaknin, R. W. Hill, Cyril Proust, C. Lupien, and L. 
Taillefer, and P. C. Canfield, Phys. Rev. Lett. 87, 237001 
(2001). 

12 E. Boaknin, M. A. Tanatar, J. Paglione, D. Hawthorn, F. 
Ronning, R. W. Hill, M. Sutherland, L. Taillefer, J. Sonier, 
S. M. Hayden, and J. W. Brill, Phys. Rev. Lett. 90, 117003 
(2003). 

13 H. Shakeripour, C. Petrovic, and Louis Taillefer, New Jour- 
nal of Physics, 11, 055065 (2009). 

14 M. Sutherland, N. Doiron-Leyraud, L.Taillefer, T. Weller, 
M. Ellerby, and S. S. Saxena, Phys. Rev. Lett. 98, 067003 
(2007) 

A. V. Sologubenko, J. Jun, S. M. Kazakov, J. Karpinski, 
and H. R. Ott, Phys. Rev. B, 66, 014504 (2002). 

16 M. A. Tanatar, J.-Ph. Reid, H. Shakeripour, X. G. Luo, 
N. Doiron-Leyraud, N. Ni, S. L. Bud'ko, P. C. Canfield, R. 
Prozorov, and L. Taillefer, Phys. Rev. Lett. 104, 067002 
(2010). 

17 J.-Ph. Reid, M. A. Tanatar, X. G. Luo, H. Shakeripour, 
N. Doiron-Leyraud, N. Ni, S. L. Bud'ko, P. C. Canfield, 



R. Prozorov, and L. Taillefer, Phys. Rev. B 82, 064501 
(2010). 

18 Y. Machida, S. Sakai, K. Izawa, H. Okuyama, and T. 
Watanabe, larXiv:1009.2432l 

19 C. Caroli, P. -G. de Gennes, and J. Matricon, Phys. Letters 
9, 307 (1964); C. Caroli and J. Matricon, Physik Konden- 
sienten Materie 3, 380 (1965); J. Bardeen, R. Kiimmel, A. 
E. Jacobs, and L. Tewordt, Phys. Rev. 187, 556 (1969). 

20 L. P. Gor'kov, Zh. Eksperim. i Teor. Fiz. 37, 833 (1959) 
[English transl.: Soviet Phys.— JETP 10, 593 (I960)]; E. 
Helfand and N. R. Werthamer, Phys. Rev. 147, 288 (1966). 

21 A.I. Larkin and Yu. N. Ovchinnikov, Zh. Exp. Teor. Fiz. 
68, 1915, (1975) [Sov. Phys. JETP 41, 960 (1976)]. 

22 J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

23 W. Belzig, F. K. Wilhelm, C. Bruder, G. Schon, and A. 
D. Zaikin, Superlattices and Microstructures, 25, 1251 
(1999). 

24 K.D. Usadel, Phys. Rev. Lett. 25, 507 (1970). 

25 V. Chandrasekhar, In: The Physics of Superconductors, 
ed. by K.H.Bennemann and J.B.Ketterson, pp. 55-1 10, 
Springer (2004). 

26 R. Watts-Tobin, L. Kramer, and W. Pesch, J. Low Temp. 
Phys. 17, 71 (1974). 

27 A. A. Golubov and M. Yu. Kupriyanov, J. Low Temp. Phys. 
70, 83 (1988). 

28 S. Dukan, T.P. Powell, and Z. Tesanovic, Phys. Rev. B 66, 
014517, (2002). 

29 X. X. Xi, Rep. Prog. Phys. 71, 116501, (2008). 

30 A.E.Koshelev and A. A. Golubov, Phys. Rev. Lett. 90, 
177002 (2003). 

31 A. Gurevich, Phys. Rev. B 67, 184515 (2003); A A. Gol- 
ubov and A.E.Koshelev, Phys. Rev. B 68, 104503 (2003) 

32 Y. X. Wang, T. Plackowski and A. Junod, Physica C 355 
179, (2001); F. Bouquet, R. A. Fisher, N. E. Phillips , D. 
G. Hinks and J. D. Jorgensen, Phys. Rev. Lett. 87 047001 
(2001). 

33 A. Shibata, M. Matsumoto, K. Izawa, Y. Matsuda, S. Lee 
and S. Tajima, Phys. Rev. B 68, 060501(R) (2003) 



